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Abstract 

By performing a fully three dimensional Hartree-Fock calculation with use of the Skyrm forces, we 
demonstrate possibility of exotic deformations violating both the reflection and the axial symmetries of 
N = Z nuclei in A ~ 60 — 80 mass region. The Y32 tetrahedral shape predicted in excited 80 Zr arises 
from a shell gap at N, Z = 40 which is enhanced for the tetrahedron deformation. Softness toward the 
Y33 triangular deformation of the oblate state in 68 Se is also predicted. 



*A talk presented by M.M. at Nuclear Structure '98, Gatlinburg, August 10-15, 1998 



1 Introduction 



Octupole deformations that violate the reflection symmetry have attracted many attentions recently in 
studies of nuclear structure JD. Although axial octupole deformation is well established in actinides and in 
neutron rich Xe and Ba region, there seems no experimental evidence of more exotic shapes that violate 
both the reflection and the axial symmetries, except for light alpha-clustering nuclei such as 12 C [Q. 

Non-axial octupole deformations have been predicted theoretically for actinides (^] and in superde- 
formed nuclei j|, ||. Instead we point out possibility of the non-axial octupole shapes in N = Z A ^ 60—80 
region of the nuclear chart || fjj . For N = Z nuclei the deformation driving shell effects which arise from 
both protons and neutrons cooperate coherently. The shell effect is indeed large in this mass region as is 
illustrated by the prolate-oblate shape coexistence in Kr and Se isotopes || and the sudden onset of large 
prolate deformation for N, Z w 38 — 40 || . Furthermore octupole instability due to the single-particle 
structure in nuclear mean- field is predicted for N, Z ~ 34 (in addition to 56, 90, 134) (k], (nj for which 
the coupling between g g / 2 and p 3 / 2 orbits are mostly responsible. 

It is found recently that the shell structures plays important roles for aggregates of metallic atoms 
(metal clusters) . In this system the mean-field for the valence electrons resembles with the nuclear mean- 
field except for its effectively zero Is potential ]l2] ]. Non-axial octupole deformatios seem to be one of 
the dominant deformations in this system as demonstrated by means of schematic models |l3|, |l4| and 
more recent mean-field calculations . In particular, the tetrahedron shape or the Y 32 deformation 

is predicted to be stable because of the shell gaps generated at N electron — 40, 70, 112, 156 [|l3], [l5| [T(|. In 
this paper, we will show that the tetrahedron shape is expected in the nucleus 80 Zr with N = Z = 40, 
being in parallel with the metal cluster system in spite of the difference in the Is potential and the pairing 
correlation present in nuclei. Furthermore, non-axial octupole deformations besides Y32 are important 
deformation modes in N = Z nuclei in A ~ 60 — 80 mass region. See also the previous publication . 



2 A fully three dimensional Hartree-Fock calculation 

For investigating non-axial reflection asymmetric deformations by means of the mean-field theory, we 
have to exclude any symmetry assumptions on nuclear deformations. Furthermore, it is important to 
use a description which allows arbitrary deformation when we search unknown shapes. From this view 
point, we adopt the Hartree-Fock method and perform a fully three dimensional calculation using the 
Cartesian mesh representation without any requirements on nuclear shapes and its symmetries. The 
Skyrm force is used as the effective nucleon-nucleon force. The pairing correlation is treated by means 
of the BCS method for the seniority force with a cutoff of the single-particle orbits, which is chosen the 
same as Ref. [|l7| [ll|. The imaginary time method is used for the iteration procedure JlSfl ■ Thus the 
basic formulation is in parallel with that of Ref. |17j , except for the symmetry treatment (they assumed 
the reflection symmetries) and calculational details. A new code has been developed independently so 
that it can also include the parity projection and full variation after projectionpcfl. The parity projected 
Hartree-Fock with Skyrm force, but without BCS treatment, has been applied to light nuclci^O] (but not 
to the present investigation). 

For the results presented below, we use the 3D mesh which are enclosed by a sphere of radius R = 13 
fm, and the mesh width of 1 fm. The parameter set of SIII [^l| is used for the Skyrm force in most 
calculations while SkM* |22| is also employed for comparison. Concerning the pairing force, we use 
G p = 16.5/(11 + Z) MeV for protons as given in Ref. 0. For neutrons, we use G n = 16.5/(11 + N) 
MeV pjj since it is natural to take the same value as protons for N = Z nuclei. In order to search not 
only the ground state but also local minimum states, we performed different runs of iterations starting 
with different initial conditions. To characterize deformation of the obtained solutions, we use the mass 
multipole moments, 

_ M$IEfrfo„,(«)|$) , 
«im= ,(m - -l,- ■ ■ ,1), (1) 
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where A is the mass number and R = 1.2A 1 / 3 fm. Here Xi m is a real basis of the spherical harmonics, 

Xl\ m \ 

where the quantization axis is chosen as the largest and smallest principal inertia axes for prolate and 
oblate solutions, respectively. We put the constraints ot\ m = 0(m = — 1, 0, 1) for the center of mass, and 
oi2m = 0(m = — 2, — 1, 1) for the principal axes. This is done by adding the constraining fields —Xi m r l Xi m 
to the mean field, where the constraining Lagrange multipliers are determined so that the constraints are 
satisfied in each iteration stepppfl. 

For the quadrupole moment, we use ordinary (/3,7) notation, i.e., a 20 — cos 7, 0:22 = ,5 sin 7, mapped 
in the /3>0,0<7<7r/3 section. To represent magnitude of the octupole deformation, we define 
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Ps = (J2 a Sm)K fom = (<4 m + <4_ m )l (m = 0,1,2,3). (3) 
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= Y l0 , 

= -J={ Y l-\m\ + Y *-\m\)> 

= ~m^Xl\m\ - Y l*m\)i ( 2 ) 



3 Shapes of ground and local minimum states 

The HF ground states and local minimum solutions for the even-even nuclei 64 Ge, 68 Se, 72 Kr, 76 Sr 
and 80 Zr are listed in Table 1. The result is essentially the same as the previous calculation || while 
the convergence of the HF iteration is taken care of more carefully in the present calculation. It is 
noted that many of these nuclei have local minimum solutions with low excitation energy, suggesting 
presence of shape coexisting excited states. The ground state shape changes from the triaxial ( 64 Ge), the 
oblate ( 68 Se, 72 Kr), to the strong prolate shape ( 76 Sr, 80 Zr) as TV, Z increases. This is consistent with the 
experimental trends || Hi- 

The reflection asymmetric solutions with non-zero octupole deformation (03 > 0) are obtained as the 
ground state of 68 Se, and as the shape coexisting excited states in 76 Sr and 80 Zr. These solutions violate 



Table 1: The ground and local minimum HF solutions obtained with the SIII force. The number in the 
first line is the excitation energy in MeV measured from the ground state solution. The second and the 
third line list the quadrupole and octupole deformations of the solutions. 





Oblate > 0, 30° < 7 < 60° 


= 


Prolate > 0, 0° < 


7 < 30° 


64 Ge 


g.s. 

/?,7 = 0.28,25° (triaxial) 
03 = 0.00 


68 Se 


g.s. 

/3,7 = 0.28,60° 

03 = 033 = 0.14 




0.25 
/3,7 = 0.25,0° 

03 = 0.00 




72 Kr 


g.s. 

/?,7 = 0.35,60° 

03 = 0.00 




1.47 
/3,7 = 0.40,0° 

03 = 0.00 




76 Sr 


2.72 
0,j = 0.14,60° 
3 = O.13,0 a2 = 0.12 






g.s. 

/3,7 = 0.50,0° 

03 = 0.00 


80 Zr 


1.56 
/?,7 = 0.19,60° 
03 = 0.00 


0.96 

0,-1= 0.00,0° 
03 = 032 - 0.23 




g.s. 

/3,7= 0.51,0° 

03 = 0.00 
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not only the reflection symmetry but also the axial one. The 68 Se ground state has finte Y33 deformation 
(/?33 = 0.14) together with oblate quadrupole deformation. The density profile shows a triangular shape, 
as seen in Fig.l. In 80 Zr, the first local minimum solution with the excitation energy of 0.96 MeV has large 
Y32 deformation (^32 = 0.23) whereas the solution has no quadrupole deformation. It has a tetrahedral 
shape as shown in Fig. 1. Although this solution violates both the reflection and the axial symmetry, 
the tetrahedral symmetry (T^) of the point group emerges as a spontaneous symmetry. Indeed the 
calculated single-particle energy spectrum shows a four-fold degeneracy which is a characteristic feature 
of the Fermion mean-field with the Td symmetry [[| (See also Fig. 3(c)). 
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Figure 1: The density contour surface at the half central density for the triangular ground state solution 
in 68 Se (left), and the tetrahedral first local minmum solution in 80 Zr (right), listed in Table 1 obtained 
with the SIII force. 



4 Non-axial octupole deformations 
Potential energy curves 

To evaluate softness of the obtained solutions toward the octupole deformations, we calculate the poten- 
tial energy curve as a function of the deformation parameters 03 m for all the independent components 
m = 0,1,2,3 of octupole deformations by means of the constraint Hartree-Fock calculation. For this 
calculation, we introduce additional nine constraints for the quadrupole deformation 0120,0:22 (equivalent 
to /3, 7) and the octupole deformations a^ m (m = —3, ..,3). The values of (3, 7 is fixed to the ones of a HF 
solution under consideration, and one of 03 m (m = 0, 1, 2, 3) is varied while other 03 m 's are fixed to zero. 
The calculated potential energy curve is plotted in Fig. 2. 

Tetrahedron shape and associated shell gap in 80 Zr 

The tetrahedral solution in 80 Zr corresponds to the minimum at 032 = 0.23 of the 032 potential energy 
curve with the spherical constraint {(3 = 0). The energy gain due to the tetrahedron deformation or 
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equivalently the energy difference between the minimum at CH32 = 0.23 relative to 0132 = is as large 
as 0.70 MeV. It is noted that the potential energy curves for a^ m {m = 0, 1,3) do not show such deep 
minimum. The reason of this can be seen in the single-particle energy diagram, Fig. 3, plotted as a 
function of the octupole deformations a^ m . Figure 3 indicates that the single-particle energy spectrum 
accompanies a large gap at N, Z — 40 which increases as increasing the tetrahedral deformation CX32 
while the other octupole deformations m = 0,1,3 do not have this feature. This is the gap specific to 
the tetrahedron deformation, and presence of such gap is due to the high symmetry (the T<j symmetry 
of the point group) (l^ . Note that the tetrahedron shell gap at particle number 40 is found also in the 
metal cluster systems 16 for which there is no Is potential. Figure 3 indicates that the tetrahedron 



shell gap emerges even under influence of the large Is term in nuclear potential. It is also noted that 
the tetrahedron minima is obtained with including the pairing. The realistic Hartree-Fock calculation 
thus supports the correspondence between nuclei and clusters suggested in Ref. |l^] on the presence of the 
tetrahedron shell gap. 




Figure 3: Neutron single-particle energies as a function of octupole deformation a^ m for m = 0,1,2,3 
((a),(b),(c),(d), respectively) calculated for 80 Zr with (3 = 0. The quadrupole constraint is fixed to (3 = 0. 
The SIII force is used. The proton spectrum is almost the same as neutrons. 
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Triangular Softness in 68 Se 



Concerning the oblate triangular solution of 68 Se, the minima is quite shallow with respect to the 033 
direction. The energy gain due to the Y33 deformation (the energy difference between the minimum and 
the solution with /?3,a33 = ) is just 0.15 MeV. This indicates that the oblate state is extremely soft 
toward the ^33 deformation rather than a rigid triangular deformation. It is noted that the potential 
energy curves for a^m = 0(m = 0, 1, 2) are not as soft as the 033 curve. 



Systematics 

In Fig. 2 and Table 1 are seen systematic trends of the non-axial octupole deformations in N = Z, A ~ 
60 — 80 nuclei. The octupole softness is enhanced at N, Z ~ 34, (representative is oblate ground state 
of 68 Se which is soft for triangular deformation), and at N, Z ~ 40 (the tetrahedral deformation). It is 
also seen that the I33 mode is the softest among the four octupole deformations for the oblate states 
whereas Y31 and Y30 is favored for the prolate states. This is explained in terms of the single-particle 
shell structures near the Fermi surface^. For the oblate states, the high f2 orbits such as [404]9/2 and 
[413]7/2 stemming from g g / 2 are located near the Fermi surface of N, Z ~ 34 and strong Y33 coupling 
with [301]3/2 and [310] 1/2 orbits emerges. On the other hand the 39/2 low orbits are situated far 
above the Fermi surface. Thus Y30 Y31 couplings are disfavored. For the prolate states Y31 is favored 
in a reversed reasoning. The Y32 deformations favored at N, Z ~ 40 is caused by the tetrahedron shell 
gap discussed above. The instability and the softness for the non-axial octupole deformations are quite 
contrasting with the axial octupole instability known in other mass regions. 



5 Discussion 



Sensitivity to force parameters 



Besides SIII, other parameter sets such as SkM* 22 are known to provide a reasonable description of 



nuclear deformations. To check sensitivity to the force parameters we performed a calculation using 
SkM*. The calculated potential energy curves associated with an oblate HF solution in 68 Se and a 
spherical HF solution in 80 Zr are shown in Fig. 4 (a) and (c). Although they do not have a minima 
at finite a^m, the potential curve shows softness for the tetrahedral 032 direction of the 80 Zr spherical 
state ((3 = 0) and for the triangular 0:33 direction of the 68 Se oblate state. This qualitative agreement 
with the SIII results can be expected since the single-particle structures which play a central role for the 
mechanism for the non- axial octupole deformations are not very different for both parameter sets. In 
Fig. 4(b) and (d), we show the results of SIII but with increased pairing force strength in order to show the 
sensitivity to the pairing correlation. Here we assumed 10% increase of G„ iP which is within a reasonable 
range reproducing the experimental pairing gap (the odd-even mass difference). From comparison with 
Fig. 2, it is seen that the pairing correlation have an effect to reduce the non-axial deformation. Note, 
however, that the essential features of the potential curve remain the same. 



Spectral signatures 

If a rigid tetrahedron deformation is realized, a characteristic pattern of excited spectra is expected. Rota- 
tional excitation of a even-even tetrahedron nucleus will accompany a sequence of levels + , 3 _ , 4 + , 6 + , 7 — , .. 
following the rotational energy relation E(I) — E(Q) = 1(1 + X)/2J in parallel with the tetrahedron 
molecule ^6|. Here the spectrum is labeled by the irreducible representation of the Td group. The rota- 
tional levels are connected by strong E3 transitions instead of E2. Viewing the potential energy curve for 
the tetrahedral solution in 80 Zr, we can expect not a rigid deformation but rather a transitional situation 
between the tetrahedron rotation and the octupole vibration around the spherical state. From the experi- 
mental systematics of the first 3~ excitation energy in the neighboring mass region |27j], a typical octupole 
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Figure 4: Potential energy curve calculated with the SkM* force (a,c), and with SIII and an increased 
pairing force G n , P = 18.15/(11 + N, Z) MeV (b,c), for the excited spherical state in 80 Zr (a,b) and the 
oblate state in 68 Se (c,d). 



vibrational energy is about a few MeV. Using this estimate and the calculated potential energy curve 
for the 0:32 direction, amplitude associated with the transitional tetrahedral excitation is estimated to be 
about a 32 (/3 3 ) ~ 0.3. This corresponds to fl(E3) - 60 W.u. on the basis of B(E3) = (3Zei? 3 /3 3 /47r) 2 , 
which is much larger than the systematic values pS| in neighboring mass region. 

For the oblate state in 68 Se, the potential energy curve suggests soft vibrational excitation for the Y33 
direction which couples with the rotational excitation. We can expect for such case that the excitation 
spectra showing a K"" = + rotational band (0 + , 2 + , 4 + , ..) associated with the ground state, and a 
K 77 = 3~ rotational band (3~,4~,5~, ..) associated with the vibrational excitation. From the calculated 
potential energy curve, the vibrational amplitude is estimated as 033(^3) ~ 0.2 — 0.3, corresponding to 
5(E3) ^20 — 40 W.u. for interband transitions. 

The nuclei under discussion are close to the proton drip line, and experimental information on the 
excited spectra is not very rich at present. Candidates of negative parity states are found in 64 Ge, which 
suggests presence of some but not very strong octupole collectivity (24). Quite recently candidate negative 
parity states are found in 68 Se, and the observed spectra looks similar to that in 64 Ge although the 
experimental information is not enough to see the octupole collectivity psf. 
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6 Conclusions 



The fully three dimensional Hatree-Fock calculation suggests that the non-axial octupole deformations 
are important in N = Z, A ~ 60 — 80 nuclei. Prominent examples are the tetrahedral deformation of 
an shape coexisting excited state in 80 Zr, and the triangular softness of the oblate state in 68 Se. The 
tetrahedron deformation causes a shell gap at N, Z = 40, which founds the microscopic origin of the 
exotic shape. 
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